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Abstract 

We consider the Josephson Effect between two spatially separated Bose- 
Einstein condensates of atoms each of which can be in two hyperfine states. 
We derive simple equations of motion for this system closely analogous to the 
Bloch equations. We also map the dynamics of the system onto those of a 



classical particle in a well. We find novel density and spin modes of oscillation 
and new stable equilibrium points of the motion. Finally we analyze the 



oscillation modes in the spin-1 (F=l) case. 



I. INTRODUCTION 

Recent experiments on the Josephson effect in 3 He have seen new phenomena that have 
not been previously observed in conventional superconducting junctions or in 4 He superleaks 
[T|. Some authors have attributed these effects to the spin of the Cooper pairs since they are 
paired in an S=l configuration ||[|. By analogy this suggests that the simple Josephson 
effect in Bose condensed alkali gases (BEC) could be qualitatively modified by the presence 
of an internal spin degree of freedom. Recently several experimental groups have succeeded 
in forming this kind of condensate in which the atoms can be in more than one hyperfine 
state. The JILA group [f| has trapped the \F = 2,rriF = 1) and \F = l,m,F = —1) states of 
87 Rb using magnetic fields. The Ketterle group at MIT || has trapped the F = 1 multiplet 
of 23 Na with optical methods. Both groups have obtained condensates with densities around 
10 14 cm -3 . An important feature of their setup is the possibility of imaging each species 
separately, thus allowing them to observe their individual motion. There have also been 
several theoretical studies of these systems |] || . 

The external Josephson effect, i.e. between two spatially separated condensates in the 

case of a single hyperfine state, has already been addressed extensively in the literature 

10| ]T6[ . The internal Josephson effect (between hyperfine states) has also been analyzed 



|T7] . In the present work we study only the external Josephson effect, that between two 
spatially separated condensates in a double well potential whose atoms have two possible 
internal states, |1) and |2). A weak link is established between the condensates by lowering 
the potential barrier that separates them. 

This paper is organized as follows: first we describe the model and the Hamiltonian of 
the system in Sec. II. After that, in Sec. Ill we derive the equations of motion and show 
that, in the so-called isotropic limit, they reduce to a form which is equivalent to the well- 
studied Bloch equations. In Sec. IV we identify the equilibrium points of the motion and 
study their dynamic stability. Finally, in Sec. V we make some considerations regarding 
the extension to the F — 1 case (when the trapping potentials for all three mp sublevels are 



identical) . 

II. MODEL SYSTEM AND HAMILTONIAN 

Let us take a condensate in a symmetric double well potential where the barrier is much 
larger than the chemical potential and therefore, in order to go from one side to the other, 
the atoms must tunnel under the barrier. Each atom has two possible hyperfine states, which 
means that the order parameter is a two component function. This setup can be achieved 
by taking a condensate in a single well and raising a potential barrier in the middle, thereby 
splitting it into two parts. Following this we may apply laser pulses to each side selectively 
in order to choose a particular superposition of internal states of the atoms on each side. 

We shall now make a four mode approximation to describe the system. Let |1,-R), 
|2, R), \1,L) and |2, L) be the four single atom states corresponding to the four modes, 
where the labels R and L refer to the right and left wells. The single atom states are, in 
principle, time-dependent and will be approximately given, in the adiabatic approximation, 
by the Gross-Pitaevskii groundstate, which in turn is determined by the number of particles 
in each single atom state. When each of the four states defined above is macroscopically 
occupied, the condensate wavefunction ^ l RL (i = 1,2) inside each of the wells can be well 
described at the Gross-Pitaevskii level. 

We shall assume that the system is always in the semiclassical regime, i.e. that the 
fluctuations around the mean values of the physical quantities are small (see below). With 
this proviso we can describe the system in terms of classical (c-number) canonically conjugate 



variables [ 19 1 . The semiclassical variables are: AiVj - one half of the difference in the number 
of atoms in each internal state, and the relative phases Atpi, which are canonically conjugate 
to each other. The condensate wavefunction can be also expressed in terms of these variables 



as ty 3 RL = \/Nj ± AA/j exp (±iA(fj/2) up to an overall phase, where the N/s are one half 
of the total number of atoms in state j. Since we are assuming in this paper that there 
is no laser coupling between the different spin states, it follows that the different N/s are 



conserved separately. Therefore the dynamics is not naturally described by using as variables 
the relative phases between internal states in each well. Although the difference in chemical 
potential between the two species in each well can be very large, this does not affect the 
Josephson dynamics. For the same reason, also the dephasing between states of different 
hyperfine spin is irrelevant. 

Furthermore, if we restrict ourselves to values of (AiVi + AA r 2)/(A f i + iV 2 ) <C f , we can 
write down an approximate Hamiltonian which is the straightforward generalization of the 
spinless case and reduces to it when one of the iVj's (and therefore the AA^'s) goes to zero: 

H = Hj + H int (f) 



where 



Hj = -u ]T y/Nf-ANfcosAipi (2) 



i=l,2 



and 

H int = l - (e 11 AN 2 1 + e 22 AiV 2 2 + 2e 12 AJV 1 AJV 2 ) (3) 

Hj is the Josephson coupling Hamiltonian and H mt is the interaction Hamiltonian (for 
a derivation of these see [jll||l5ll ). The interaction term conserves the total number of atoms 
in each hyperfine state separately and we are here ignoring loss processes that occur in the 
real system, ujq is the Josephson tunneling energy which we take to be the same for both 
hyperfine states for simplicity; and e^- are the effective interaction coefficients. We assume 
that ujq is independent of AiVj and Ay?j, although we do expect some dependence for large 
values of AiVj. We shall also assume that all e's are positive because we want to avoid two 
possible complications: the collapse of the gas and a possible accumulation of atoms on one 
side of the junction. 

In the spinless case we can identify four different dynamical regimes as the parameters 
ujq and e are varied. By analyzing the topology of the phase space (Fig. la-b-c), we find 
two sharp transitions and one crossover. Starting from the weakly interacting limit (the 



so-called Rabi regime - Fig. la), the transition to the intermediate regime is marked by 
the appearance of closed orbits oscillating around nonzero values of AiV and centered at 
Aip = it (Fig. lb). The second occurs when open orbits appear where Acp extends over all 
values (Fig. lc). This is a transition from the intermediate to the Josephson regimes. These 
transitions happen at the values Uq = eN and 2uo = eN respectively for our model. Finally 
when ujqN ~ e there is a crossover to the Fock regime where quantum phase fluctuations 
cannot be neglected. Deep in the Rabi regime (when ujq ^> eN) the tunneling energy 
dominates. In the Josephson regime both Hj and H int are important (u <C eN <C u>oN 2 ) 
and finally, in the Fock regime, H[ nt dominates (ujqN <C e) [ |l9|j . 

u can be varied anywhere from to 100 s _1 . On the other hand, e can go from 0.01 
to 0.1 s _1 |L5j and, under current experimental conditions, N is usually between 10 4 and 



10 7 . With this range of parameters the Fock and Josephson regimes are easily accessible 
whereas the Rabi regime is more difficult to achieve. It is important that the frequency of 
any oscillation between the wells be smaller than the lowest intrawell excitation frequencies 
so that, during the motion, these degrees of freedom are not excited. In practice this means 
that both ujq and ^/uj^eN have to be smaller than the frequency of the lowest intrawell 
collective mode (as will become clear below). 

To be consistent with the semiclassical description we require that the standard de- 
viations of the quantum operators A(p and AiV satisfy the conditions cr(Aip) <C 1 and 
a (AN) <C N during the motion of the system. Generally speaking, the experimental setups 
will be such that cr(AN) ~ o~ l (A0) ~ (e/cuoN + l/N 2 )' 1 ^. The second inequality is always 
satisfied for positive e's. The first is satisfied only in the Josephson and Rabi regimes, to 
which we shall restrict our analysis from now on. 

In order to justify the effective Hamiltonian the two components must be miscible, in 
other words, there can be no component separation. If this were not the case we might not 
have the same tunneling matrix element o>o for both species and the mean field interaction 
energy would not have the form that we assume. This means that a condition must be 



imposed on the interaction parameters, namely that ene22 > ef 2 |18 



Finally, the experimental observations can be made by measuring the density (and there- 
fore AA^ 12 ) in the usual way, either destructively or by phase-contrast imaging. As men- 
tioned before, an important point is that these methods allow us to determine experimentally 
the behaviour of each hyperfine species separately. 

III. DYNAMICS AND BLOCH EQUATIONS 

The equations of motion are: 



AN, = - ^— = -u y/N? - AN? sin A Vl (4) 

fjlT A /V 

A Vl = — = u , _ ' cos A^ + £ eyANj (5) 

\ 1 % J 



In this section we will rewrite the equations of motion in terms of new variables in order 
to provide some insight into the dynamics of the system. 

A. Isotropic case 

We first consider the isotropic case where en = e\ 2 = £22- We shall be working with the 
quantity 

e = -(e u + e 22 + 2e 12 ) (6) 

which is in fact equal to any of the e's in the isotropic case. This definition however, 
will be useful ahead, when we deal with the anisotropic situation. The equality of the 
interaction parameters in the isotropic case seems to violate the miscibility condition that 
£11^22 > €12 ■ However this condition does not take into account the kinetic energy which 
favours miscibility. Therefore isotropy does not pose any such problems. 

We notice that the Hamiltonian (|1|) is invariant under arbitrary spin rotations and more 
generally of £77(2) transformations applied simultaneously to the spins in both wells. That 
is, if we transform the two-component spinors ^ l and ^> r with the same unitary operator 
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then the dynamics should remain unchanged. This suggests that we re-express the equations 
of motion in terms of quantities which are invariant under such transformations. This con- 
clusion of course depends on the isotropy of the interaction Hamiltonian H mt . We therefore 
define the following dot products of spinors: 



AiV, EE 



KTr 12 |,Tr 12 

V^R\ — \*L\ 



J2AN t 



a. 



V* L V R - y* R y L 



2i 



J^^Nf-ANismAipi 



P+ = 



^H 



**R*L 



J2 \/N?-ANi cos Aw 



(7) 
(8) 
(9) 



The subscript (+) will be used to distinguish this set of variables from another one with 

subscript (— ) to be defined below. Using the equations of motion for AiVj and Aipi, we 

obtain 

/AN + \ / -u \ /AN + \ 



(i+ = w +eA7V + 

\ j3 + J \ -eAiV + 

If we now define the three-component vectors i% 



(10) 
(AiV + ,a + ,/3+) and B(t) = 



(— eAN + , 0, uq) then we can rewrite the equations of motion succinctly as 



B(£)xr H 



(11) 



Note though that B and r + are not independent since they are both functions of AiV + . 
Straightforward manipulation of Eq. ([TOD, or directly of Eqs. (H) and (Bl), leads to 



A'N + = - (ul - 6/7(0)) AN + - "-ANl 



(12) 



where H(0) = -u P+(0) + f AiV + (0) 2 , (N = J2>N t ) 

This equation is quite general since it is valid not only for two hyperfine states but for 
any number of them, as long as they interact only through a AN+ = (J2i AiVj) 2 term. In 
particular it also applies to a single state (i.e. spinless) system. It is formally identical to 
the equation of motion of a particle with unit mass in the quadratic-plus-quartic effective 
potential 



V eff (AN + ) = \ (ujI - e /7(0)) ANl + ^ANl (13) 

with effective total energy 

E eff = V ef ,(AN+) + l -AN + 2 (14) 

An important point to notice is that E e ff and V e ff cannot be chosen independently 
since they both depend on the initial conditions. The variation of H(0) and of a + (0) (since 
AN = — uoa+) allows us to find three different types of motion (Fig. 2a-b-c). In the first 
type the coefficient of the quadratic term is positive and AN + oscillates around zero, which 
is the minimum of V e ff (Fig. 2a). In the spinless case this corresponds to either oscillations 
around the origin (Fig. la-b-c) or to small oscillations around the 7r-state (Fig. la). The 
second case occurs when the coefficient is negative and E e ff is positive, which also leads to 
oscillations of AN + around zero although that point is no longer a minimum of V e ff (Fig. 
2b). It corresponds to large oscillations around the origin or the 7r-state (Fig. lb-c). The 
third one corresponds to both E e ff and the coefficient being negative (Fig. 2c) and leads 
to self-trapped behaviour (oscillations around AN + ^ - Fig. lb-c). For the spinless case 
there is a well-known analogy with a momentum-shortened pendulum in a gravitational field 
whose behaviour is also fully reproduced by this particle-in-a-well model. Since the analysis 



of the spinless junction has already been carried out in Ref. [Oj we shall not continue it 
here and shall proceed to the two hyperfine state case. 

Specifying the dynamics of AN + does not describe the motion completely. For example, 
even in the spinless case it is known that the third regime includes two different behaviours 
of the relative phases, the so-called "running" and "oscillating" phases. For a description 
of these as well as 7r-states and the momentum-shortened pendulum analogy see e.g. Ref. 



II]. To further understand the dynamics of the two hyperfine state Josephson effect we 



introduce the additional variables 

AAL = AiVi - AiV 2 (15) 



a_ = JNl - ANl sin A Vl - J N% - AiVf sin Ay 2 (16) 



/3_ = y/Nf-ANfcosAipi - ^Af -AAfcosA^ 



(17) 



and their equations of motion are: 
/AAL\ /0 

w 






\ /AM\ 

+eAiV 4 



/ 






(18) 



-a; 



V -eAiV + 

Since the matrix is the same as in Eq. ( |TUD we define the three-component vector 
r = (AAL, a_,/3_) and rewrite the equations of motion as 



r = B(i)xr_ 



(19) 



Now, however, B and r_ are independent and therefore these equations are formally 
identical to the Bloch equations, familiar from the context of NMR and quantum optics. 
Notice that in going from the original four variables to six we are enlarging the configuration 
space, which means that not all points described by the new set of variables are physically 
allowed. Therefore care must taken in choosing the initial conditions of the motion. 

We can obtain some physical insight into the variables AN± by noting that AiV + is 
one half of the difference in total number between the right and left wells and, in the limit 
N\ = N2, AAL is one half of the difference in spin between the wells (|^/^| 2 — |\l/|.| 2 )/2 — 



(I* 



1 12 



ty L \ )/2. This means that the former describes the density mode whereas the latter, 



in that limit, describes the spin mode. 

We can now analyze Eqs. (|TT|) and (|19"D in a few limiting cases to gain some insight 
into the behaviour of the system. It is possible to have small oscillations in AA" + and no 
motion in AAL or vice-versa. We consider two cases: the Rabi limit where eA^ ^C ujq and 
the Josephson regime, where eA^ ^> uj . In the Rabi case, neglecting higher order terms, 
the frequency of oscillation of AA^ + (and therefore of the density mode) can be calculated 



from Eq. ([13]) to be Juq + 0^/?+- Also, using Eq. (jI9|), we can neglect the component of 
the B field along the AAL axis so that r_ (the spin mode) rotates around the /?_ axis with 
frequency ujq. In the Josephson case we consider two types of situations - small oscillations 



of AiV + around zero and around nonzero values. For zero values we can get density and 



spin modes with frequencies J Uq + u^eN and uj respectively. For nonzero values (i.e. when 
AN + is "self-trapped" around a value AJV?), AN + oscillates with frequency e | A A^. | and 



AiV_ with frequency Jul + e 2 (Ai\^) 2 . 

B. Anisotropic case 

As we would expect, the equations in this case become much more complicated. However, 
for the sake of completeness, we include them here. We find that the equations for r + and 
r_ become coupled: 

r+ = B 1 (t)xr + + B 2 (t)xr_ (20) 

r =B 1 (£)xr_ + B 2 (£)xr+ (21) 

where Bi = (-eAA+ -e fl AJV_, 0,u ), B 2 = (-e^AAL -e B AA+, 0, 0), e = J(en + e 22 + 
2e i2 ) as before, e A = |(eii + e 22 - 2e 12 ) and e B = \{e xl - e 22 ). 

IV. DISCUSSION OF THE EQUILIBRIUM POINTS OF THE MOTION 

In this section we study the existence and stability of the equilibrium points of the 
motion. The main results for the isotropic case are summarized in Table 1. 

At an equilibrium point, ANi i2 = A<^i i2 = 0. Using Eq. (f|) this implies that the phases 
Aifi and Acp 2 are either zero or it. From Eq. (|5]) we get 

AiV? = -AiV 2 ° ( 2* + , ^ ) (22) 

V 612 6 12V /iV 2 2 -(AA 2 °)2C 2 y 



AiV 2 ° = -AN? — + ; ^° ) (23) 

> ei2 JN?-(AN?yCi 



where AN® and Aip® are the coordinates of the equilibrium point. We have defined 
£i = cos A(/? t j ) and ( 2 — cos Aip® to abbreviate the formulae. 
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These equations define two functions, AN® (AN®) and AN® (AN®), which we can plot 
on the (AN®, AN®) plane (Fig.3). Since both d and (2 can be 1 or —1 (corresponding to 
Aip® = or 7r) we have three distinct cases. For all of them the trivial point AN® = AN® = 
is always a solution. Therefore in all cases we have at least one solution. 

Case 1: & = C2 = 1 (A<^ = Ay? 2 = 0) 

As is clear from Fig. 3a-b the condition for the existence of three solutions imposes a 
condition on the slopes of the curves at the origin which leads to 

Case 2: Ci = -1, C2 = 1 (Api = tt, A^ 2 = 0) 

For three solutions to exist (Fig. 3c-d) we require this time that 

/en _ u \ /£22 ^0 \ 1 (25) 

U12 € 12 nJ U12 e 12 N 2 J { ' 

Case 3: Ci = C2 = -1 (A^i = Atp 2 = tt) 

This time we may have one (Fig. 3e), three (Fig. 3f) or five (Fig. 3g) solutions. To have 

three we need that 

'e 22 ^0 \ / e n ^o 



-.. , < ! ( 26 ) 

,ei2 ei2N 2 J Vei2 £12-^1/ 

If this condition is not met then we will have one or five solutions, depending on whether 
the factors on the left hand side are both negative or both positive respectively. 

In the rest of this section we will analyze the behaviour of the system close to the various 
equilibrium points. Notice that the global groundstate is the trivial solution AAT 1)2 = 
and Aipip = 0. All the other equilibrium points are thermodynamically unstable although 
possibly dynamically stable p3f . 

A. Isotropic case 

If all €i/s are equal then some of the conditions in the previous subsection cannot be 
satisfied. For Case 1, condition fl2~3| ) cannot be satisfied and therefore only the equilibrium 
point AiVi = AA 2 = is allowed. 

11 



In Case 2 both single and triple solutions are allowed. Condition (|25|) for the existence 
of three equilibrium points becomes Aq — N 2 — uj /e > 0. 

In Case 3 the conditions for the existence of one or three equilibrium points can be 
satisfied. The condition for three is Aq + N 2 — wo/e > 0. However we cannot have five 
equilibrium points since, if both terms are positive in condition ( PE|) and each of them is 
smaller than one, their product will also be smaller than one. 

To study the behaviour in the neighbourhood of an equilibrium point we shall work with 
the second order differential equations for AAq and AA^. To obtain these we differentiate 
Eq. (HD with respect to time and eliminate Atpi^ and AAq^ using Eqs. (|4]||). We now 
introduce the variables 61,62 defined by 

AAq = AA^ + 5i (27) 

AN 2 = AiV 2 + 6 2 (28) 

The linearized equations of motion for the isotropic case are: 

^U-fiM^ (29) 

62) w 

where 



n 2 =U + e 2 (AN° + ) 2 Id + ^oe " \ (30) 



Case 1: Ci = C2 = 1 

As mentioned before, the only stable point is at AN® = AN® = 0. In the basis (61,62) 
we find the modes (Aq, N2) and (1,-1) (note that the matrix of the linearized equations 
of motion is not hermitian and therefore the two eigenvectors are not guaranteed to be 
orthogonal even if the corresponding frequencies are different). The first corresponds to a 



'density' mode with frequency Jcu 2 + ujqeN and the second to a 'spin' mode with frequency 
ojq. Note though that, even in the density mode, the total spin on each side of the junction 
changes as a function of time (unless Aq = A^ 2 ). 

12 



Case 2: Ci = -1, C2 = 1 

Near AN® = AN® = we proceed as above and find the eigenfrequencies 



uJq + ujqe(N2 — Ni) and cu with corresponding eigenvectors (Aq, —N 2 ) and (1, —1). The 
system is dynamically stable as long as the frequencies are real, which leads to the condi- 
tion Aq — N 2 — uj /e < 0. Since N% and N 2 are easy to change experimentally, this state 
can always be made stable regardless of the values of ujq and e. It is therefore much eas- 
ier to obtain a ir state this way than in the spinless case. However there is an additional 
complication: if Aq = N 2 the eigenvectors become parallel and, since the representation of 
arbitrary vectors in an almost colinear basis can involve very large amplitudes (especially for 
the vectors perpendicular to the basis vectors), the amplitude of the oscillations for initial 
displacements in the in- phase direction will tend to diverge as (A 2 — Ai)/A — > after some 
time. Experimentally it is not difficult to avoid this pitfall. If AN® 7^ and AiV^ 7^ the 



frequencies become e|AA^| for the density mode and Juq + e 2 (AiV|) 2 for the spin mode. 
Both frequencies are real and therefore this equilibrium point is stable. 
Case 3: Ci = C2 = _ 1 



For AN® = AN® = we find the eigenfrequencies Juq — uoeN and uq with corre- 
sponding eigenvectors (Ni,N 2 ) and (1,-1). The system is dynamically stable as long as 
iV < u>o/e. If AN®, AN® 7^ then the frequencies of the two eigenmodes are given by the 



same expressions as those for AiV-f, AN® ^ in Case 2 (e| AN®_ \ and Jcoq + e 2 (AiV^) 2 ) and 
therefore they are both stable as long as the points exist. 

Note that all the frequencies found in Sec. Ill A are in agreement with those derived 
here by studying the small oscillation behaviour directly from Eqs. ([|,|5]). 

B. Anisotropic case 

In the anisotropic case we can still use Eq. (E9|) but with O 2 given by 
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ul + (enAiV? + ei 2 AiV 2 ) 2 

ujI + (e 22 AiV 2 + eisAATO) 2 



e n ^N? - (AiWG W^i - (AiV?)V , 

+ ^o , , (31) 

e 12 JNi - (AiV 2 °) 2 C 2 e 22 ^iV 2 2 -(AiV0)2C 2i 



Case 1: Ci = C2 = 1 

When AiV° = AiV 2 = the eigenvalues are 



2 2 ^ e n N 1 Ci + e 22 N 2 C 2 ^ (gniViCi ~ e 22 iV 2 C 2 ) 2 , 2 , r , r . . \ ,, . 

w = ^ + ^0 2 ± V 4 + ^2^1 ^2(1(2 (32) 

with £1 = £2 = 1. For simplicity we shall address only the nearly isotropic case. It 
is experimentally relevant since this is the case for 23 Na where the experimental values for 
the en, e 22 and e\ 2 are similar. To do this we use the variables e, €a,b defined in section 
III B since e^ and e# quantify the degree of anisotropy. We therefore treat them as small 
parameters. Expanding the square root in Eq. (5^) and keeping terms to first order in those 
variables we obtain the two eigenvalues 

u 2 = u 2 Q + w e(iViCi + iV 2 C 2 ) + uj^a ^ 1 '^ + 2u e B (N 1 ( 1 - N 2 ( 2 ) (33) 

iViC,i + 1\ 2 L, 2 

U =U} ° +4UJ0€A N 1 t 1 + N 2 t 2 (34) 

We have assumed that N\ + N 2 ~ iVi — iV 2 ~ iV. Since e^ > (which is implied 
by the miscibility condition) both modes are stable (e# <C e). The instability that would 
arise at sufficiently large and negative values of e^ has the same origin as the immiscibility 
condition. However we do not consider this region in this paper since immiscibility would 
have severe consequences (see end of section II). It is easy to see from Eq. (p4|) that the case 
AN®, AiV 2 7^ is also ruled out due to the miscibility condition. 

Case 2: Ci = -1,(2 = 1 
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At the origin AN® = AN® = the frequencies are again given by Eqs. (|33|) and fl34|) 
but with d = — 1, C2 = 1- Let us divide the region into two parts: Ni — N 2 > and 
Ni — N2 < 0. In the first region the motion is always unstable: for large values of N\ — N2 
the first frequency is imaginary and, for small values, the resonance of section IV A will tend 
to destabilize the equilibrium point. In the second region, if Ni — N 2 > — 4eANiN 2 /uj then 
again it is unstable. Otherwise it is stable (provided it is outside the region of resonance). For 
AN®, AN® j^z 0, the conditions for stability become rather complex and offer little insight. 
However, on physical grounds, by suitably choosing the parameters, any of the equilibrium 
points can be made stable. 

Case 3: Ci = C2 = — 1 

The frequencies at the point AN® = AN® = are those given by Eqs. (|33"D and fl34|) 
with £1 = — 1, £2 = — 1- As in the isotropic case, the point AN® = AN® = is unstable for 
typical experimental conditions (see below), namely, when eiV > uq. The corrections to this 
criterion are of order (e^, egj/e. Finally, when AN®, AN® 7^ 0, as in the preceding case, the 
stability can generally be achieved for all equilibrium points for appropriate values of the 
parameters barring immiscibility problems. 

C. Experimental considerations 

The typical frequencies of small oscillations can be calculated using the following param- 
eters: N ~ 10 6 atoms, e ~ 0.01 s _1 , 6a,b ~ 10~ 4 s _1 and uj ~ 10 s~\ For these values most 
of the frequencies lie between 10 s _1 and 100 s~ l , whenever stable oscillations exist. 

For a general initial state near the trivial equilibrium point, which is Atpi^ = and 
ANip = 0, the oscillations in AAq^ and Atpi^ will be a superposition of both density and 
spin modes. For the typical parameters that we are using, the frequency of the density 
mode is one order of magnitude larger than that of the spin mode and, therefore, it should 
be simple to distinguish between them experimentally. It should also be possible to prepare 
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an initial state in which only one of the two modes is significantly excited. 

Although the ^-states are unstable in the spinless case for typical parameters, we have 
shown that they can be stabilized in spinor condensates. To prepare them experimentally 
we must have N\ < N 2 as explained above, where the 7r-phase difference is in species 1. 
A frequency measurement of the density mode could be used to detect that in fact a 7r 
phase exists in species 1. Alternatively, one could observe the destabilization of the state 
suddenly appearing in the form of density oscillations due to the reduction of N 2 . Finally, 
a third possibility would be the direct imaging of the interference pattern between the left 
and right condensates of species 1 during its expansion, after the trapping potentials have 
been switched off. 

V. F-l SPIN JOSEPHSON EFFECT 

In this section we look at the mean-field groundstate and Hamiltonian of a Josephson 
junction containing atoms with F — 1 total spin. The groundstate of a single spinor conden- 
sate has been analyzed in the literature [§-|9| . Here we extend the analysis to the case where 
the condensate is comprised of two spatially separated parts linked by a weak junction. 
Some of the results in this section are similar to those of Ref. || and are related to the bulk 
excitation spectrum of spin-1 condensates. As in the previous sections we assume that the 
trapping potentials are identical for all three hyperfine states (which can be achieved using 
optical dipole traps). Under these conditions it is known that all three hyperfine states of 
the multiplet are miscible. 

One might try to proceed as in section III by deriving a set of equations for invariant 
quantities such as AN+, a + , [3+ and so on. However it turns out that while this is possible, it 
does not lead to simple equations of motion as in the two internal state system and therefore 
this approach does not seem to provide a clear insight into the dynamics. 

We shall now study the small oscillations around some of the equilibrium points in both 
the ferro- and anti-ferromagnetic cases. 
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H = Hj + H h 



(35) 



where 



and 



Hj = — ^ J2 a \,L a hR + h - c - 



H m t — 2^ —(a li a li ai^ai y i + a 0i a jaojao ) i + a_ li a_ li a_i y ia_ij 

i=R,L 4 



(36) 



.t J 



J J 



t J 



+2 a l,i a 0,i a 0,i a l,i + 2a ,i a -l,i a -l,i a 0,i + ^ a l,i a -l,i a ~l,i a l,i) 
+ ~r( a l,i a l,i a l,* a l,* + a -l,i a -l,i a -l,i a -l,i + ^ a l,i a O,i a O,i a l,i 

-\-2a,Q fl—i jQ— i, i a o,i H — 2a^ jO_i jQ— i,jQi,i + ^ a i i a — i i a o,j a o,; 

+ 2ao,i a O,i a l,* a -l,i ~ a lA a l,i — a -l,i a -l,i) 



(37) 



For e.g. 23 Na, e 2 is approximately 2% of e . 

We now derive the equations of motion and, in the mean-field approximation, since we 
are assuming a macroscopic occupation, we linearize by keeping only terms at least of order 
N in H int . 

For the ferromagnetic case, if we assume that only the mp = 1 has macroscopic occupa- 
tion, we obtain the equations: 

fS(f> hL \ ( S(f>x, R \ fe + e 2 \ (25<P hL + 5<P\ L \ 



d 



i — 
dt 



0(Po,l 



Uq 

' 2 



090,R 

\S(f>-i, R J 



N 





V o 



84>o,i 



(38) 



/ 



eo + e 2 

e -e 2 J V &4>-x,l 

and a similar set for <50j,r. 

Solving them gives the following results: for the groundstate (the relative phase between 

i t n and <pi t L equal to zero), we have the following three modes: a density mode with 



frequency Jul + u> (e + e 2 )iV, a spin mode with frequency uj and a quadrupole mode with 
frequency u + |e 2 |iV. 



For the 7r-state we find the same modes with frequencies Jujq — cj (eo + e 2)^, — ^o and 
— u + \e 2 \N. The density mode can clearly become unstable for u < (e + e 2 )iV, which is 
the case for the typical parameters quoted in the previous section. 
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For the anti-ferromagnetic case, if we assume that only the mp = state is macroscopi- 
cally occupied, we obtain the equations: 



.d 

% ~dt 



I 8<f>l,L \ 



I 



UJO 

' 2 



-i,rJ 



N 



/e + e 2 \ / 

e 
V eo + W V 



h,L 



\ 






{W-iA 



+ 



Neo 







\L 



(39) 



and a similar set for 50,^. 

Solving them gives the following results: for the groundstate (the relative phase between 
t>o y R and 0o,l equal to zero), we have the following three modes: a density mode with 



frequency Juq + uoeoN and two degenerate spin modes with frequency Juq + U062N 



For the 7r-state we find the same modes with frequencies JuJq — u eoN and J Uq — u €2N. 
The density mode becomes unstable for uoq < €qN and the spin modes become unstable for 
ujq < e 2 N. With the parameters that we are using at the very least the density mode is 
unstable. 



VI. CONCLUSIONS 



In this paper we considered the Josephson junction between two spatially separated 
condensates with a hyperfine degree of freedom. We have derived a set of simple equations 
which, in the isotropic limit, are formally identical to the Bloch equations and which provide 
insight into the dynamics of the two hyperfine state condensate in a double-well setup. We 
find a partial mapping to the simple problem of a particle in a ±x 2 + x A type potential which 
becomes a complete mapping in the spinless case. We have also demonstrated the existence 
in this system of new density and spin oscillation modes. In particular we have found 7r-states 
that are stable under experimentally accessible conditions due to the interactions between 
the two species. Finally we analyzed the spin-1 case in the same geometry both for the ferro- 
and anti-ferromagnetic cases and found the low-lying oscillation modes. Our results indicate 
a wide range of phenomenology for Josephson oscillations when the superfluid has a spin 



degree of freedom. Future possible directions of research might include tunneling between 



fragmented states |H|] and more general solutions of the Bloch equations. 
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FIGURES 

FIG. 1. Orbits in phase space of the spinless Josephson effect: a) Rabi regime ujq/eN = 
1.2; b) Intermediate regime uo/eN = 0.6; c) Josephson regime ujq/eN = 0.4. 

FIG. 2. V e ff(AN + ) for different initial conditions: a) when the coefficient of the quadratic 
term in Eq. (0), namely \uJq — eH(0), is positive; b) when the coefficient is negative and 
E e ff > 0; c) when both the coefficient and E e ff are negative. The horizontal line corresponds 
to E eff . 

FIG. 3. Equations (|22]) , (|23|) plotted on the (AN\, AN 2 ) plane. The intersection of the 
two curves is the graphical solution for the equilibrium points for A(pi = Atp 2 = (a, b), 
A(fx = Ti and Aif2 = (c, d) and Atpi = Atp 2 = n (e, f, g). 

TABLE I. Equilibrium points and oscillation modes around them in the isotropic case. 
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Aip lt Aip2 


AN 1 ,AN 2 


Existence condition 


Type of mode 


Frequency 


Stability condition 


0,0 


0,0 


always exists 


density 
spin 




always stable 


Wo 


0,7T 


0,0 


always exists 


mixed 
spin 




N x -N 2 - uj /e < 
always stable 


V«g + woc(iV2-JVi) 


^o 


N 1 -N 2 - Lu /e > 


mixed 
spin 


e\ANl\ 


always stable 


^l + e^ANlf 


7T,7r 


0,0 


always exists 


density 
spin 




Ni + N 2 - uj /e < 
always stable 


^JujI - uj a eN 

LUq 


^o 


Ni + N 2 - Lu /e > 


density 
spin 


e\AN°_\ 


always stable 


^l + e^ANlf 
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